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Cooperative  Power  Scheduling  for  a  Network  of 

MIMO  Links 

Xiang  Dong,  Yue  Rong,  Member,  IEEE,  and  Yingbo  Hua,  Eellow,  IEEE 


Abstract — A  cooperative  power  scheduling  algorithm  devel¬ 
oped  hy  Wang,  Krunz  and  Cui  is  extended  for  an  ad  hoc 
network  of  MIMO  links.  This  algorithm,  referred  to  as  price- 
hased  iterative  water  filling  (PIWF)  algorithm,  is  a  distrihnted 
algorithm  by  which  each  link  computes  its  power  scheduling 
through  an  iterative  and  cooperative  process.  The  cooperation 
among  all  links  is  achieved  by  adaptive  price  factors  applied  by 
each  link.  Compared  to  a  centralized  power  scheduling  algorithm, 
the  PIWF  algorithm  is  much  more  efficient  in  computation 
although  not  as  efficient  in  network  throughput.  Compared  to  a 
non-cooperative  connter-part  by  Demirkol  and  Ingram  where  all 
price  factors  are  zero,  the  PIWF  algorithm  reqnires  additional  in- 
network  computation  but  is  more  efficient  in  network  thronghpnt. 

Index  Terms — Ad  hoc  network  of  MIMO  links,  iterative  water 
filling,  power  scheduling,  network  throughput,  cooperative  game. 


1.  Introduction 

For  scattering-rich  environment  such  as  the  case  where  the 
transmitters  and/or  receivers  are  close  to  ground,  a  MIMO 
(multi-input  multi-output)  wireless  link  is  known  to  have  a 
higher  capacity  than  a  traditional  SISO  (single-input  single¬ 
output)  wireless  link.  In  recent  years,  there  has  been  a  strong 
interest  in  developing  basic  theories  for  networking  of  MIMO 
links.  One  of  the  important  issues  for  a  network  of  MIMO 
links  is  power  scheduling,  e.g.,  see  [1],  [2],  [3],  [4]  and  the 
references  therein.  In  [1],  a  distributed  iterative  water  filling 
algorithm  was  proposed  where  each  link  tries  to  maximize  its 
own  capacity  in  a  non-cooperative  way.  This  algorithm  is  also 
similar  to  the  closed  loop  multiple-transmitter  and  multiple- 
receiver  algorithm  shown  in  [2].  In  [3],  a  centralized  algorithm 
was  developed  to  maximize  the  network  throughput  where 
the  temporal  freedom  in  time  or  frequency  was  however  not 
considered.  In  [4],  both  the  spatial  and  temporal  freedoms  were 
exploited  in  a  centralized  scheduling  algorithm. 

On  the  other  hand,  there  have  been  many  research  activities 
on  power  scheduling  for  networks  of  SISO  links.  Recent 
examples  include  [5],  [6]  and  [7],  where  several  different 
utility  functions  are  considered.  The  work  [7]  by  Wang,  Krunz 
and  Cui  is  particularly  relevant  to  the  studies  in  [1],  [2],  [3], 
[4].  In  [7],  the  authors  developed  a  distributed  algorithm, 
which  we  refer  to  as  the  PIWF  algorithm,  to  maximize  a 
priced-based  utility  function  for  an  ad  hoc  network  of  SISO 
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of  Electrical  Engineering,  University  of  California,  Riverside,  CA,  92521. 
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Department  of  Electrical  and  Computer  Engineering,  Curtin  University  of 
Technology,  Bentley,  WA  6102,  Australia.  E-mail:  y.rong@curtin.edu.au.  This 
work  was  supported  in  part  by  the  U.  S.  Army  Research  Laboratory  under 
the  Collaborative  Technology  Alliance  Program. 


links.  The  price  factors  are  chosen  in  such  a  way  that  they 
are  simple,  adaptive  and  useful  for  an  improved  network 
throughput.  In  this  letter,  we  present  an  extension  of  the  PIWF 
algorithm  for  a  network  of  MIMO  links  and  compare  its 
performance  with  [1]  and  [4].  We  will  also  refer  to  the  MIMO 
version  of  the  PIWF  algorithm  simply  as  the  PIWF  algorithm. 

Our  study  shows  that  the  PIWF  algorithm  is  indeed  more 
efficient  in  network  throughput  than  the  algorithm  in  [1], 
referred  to  as  DI  algorithm,  which  is  a  special  case  of  the 
PIWF  algorithm  with  the  price-factors  set  to  zero.  Compared 
to  the  centralized  algorithm  in  [4],  the  PIWF  algorithm  is  not 
as  efficient  in  network  throughput  but  is  much  more  efficient 
in  computation.  We  will  also  evaluate  the  network  through¬ 
put  of  the  PIWF  algorithm  using  bits-meter  per  second  per 
Hertz  per  node  (bits-meter/s/Hz/node),  which  is  a  fundamental 
throughput  unit  introduced  in  [8]  and  also  used  recently  in 
[9]  and  [10]  for  large  ad  hoc  networks.  We  will  show  that 
for  a  medium  sized  network  (64  nodes),  the  throughput  of 
the  PIWF  algorithm  can  exceed  substantially  those  in  [9] 
and  [10].  However,  the  complexity  of  the  PIWF  algorithm 
increases  as  the  network  size  increases.  For  this  reason,  the 
PIWF  algorithm  is  not  suitable  for  large  ad  hoc  networks. 

II.  Network  Model  and  Problem  Formulation 

We  assume  that  there  are  N  active  MIMO  links  that  need  to 
be  scheduled  to  share  K  frequency  bands  within  a  time  win¬ 
dow  of  interest.  Each  MIMO  link  has  Nt  transmitting  antennas 
for  the  transmitter  and  receiving  antennas  for  the  receiver. 
The  channel  matrix  between  the  transmitter  of  link  j  and  the 
receiver  of  link  i  in  band  k  is  denoted  by  Hj^i{k)  € 

The  waveform  transmitted  from  the  source  of  link  i  in  band  k 
is  Xi{k)  G  The  waveform  received  by  the  destination 

of  link  i  in  band  k  is  yi{k)  G  The  relationship  between 

these  waveforms  is  given  by 

N 

yi{k)  =  Hi^t{k)xi{k)  +  ^  Hj^i{k)xj{k)  +  Wi{k)  (1) 

i=i  jY* 

where  Wi{k)  is  the  noise  vector  which  is  assumed  to  be  white 
Gaussian  with  zero  mean  and  the  covariance  matrix  NqI.  Here, 
iVo  is  a  scalar  and  /  is  an  identity  matrix. 

The  channel  capacity  (in  bits/s/Hz)  of  link  i  in  band  k  is 
known  [12]  to  be 

r,{k)  =  log2  \I  +  P,{k)G^{k)\  (2) 

where  Pi{k)  =  E  E  is  the  statistical  expecta- 
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tion,  ^  is  the  Hermitian  transpose,  and 

G,(fc)  =  H^,{k)M-\k)H,^,{k)  (3) 

N 

M,{k)  =  NoI+  Hj4k)P,{k)H^,{k)  (4) 

i=l 

Furthermore,  if  Gi{k)  is  given  and  has  the  eigen¬ 
value  decomposition  Gi{k)  =  Qi{k)Ai{k)Q^ {k)  where 
Qi{k)  is  the  matrix  of  its  eigenvectors  and  Ai{k)  = 
diag  x\^\k)  X^i^\k)  ■  ■  ■  x[^^\k)  is  the  matrix  of  its  eigen¬ 
values,  then  the  optimal  structure  of  Pi{k)  to  maximize  the 
rate  of  link  i  is  known  [12]  to  be  Pi{k)  =  Qi{k)Ili{k)Qf  (k) 
where  ni(fc)  =  diag  7r-^^(fc)  {k)  ■■■  'n[^*\k)  . 

Given  the  optimal  Pi{k),  we  can  write  ri{k)  = 
E^:i^og,[i  +  xf\k)4\k)). 

The  main  interest  here  is  to  develop  a  distributed  algorithm 
for  computing  T:f\k)  for  all  i,k,l  so  that  the  network 
throughput  is  maximized  (at  least  locally)  subject  to  a  power 
constraint.  We  will  assume  that  there  is  a  fixed  value  for  the 
power  of  each  link,  i.e.,  we  have  such  that 
K  K  Nt 

—  <p,  (5) 

k=l  fc=l  1=1 

Although  omitted  in  this  letter,  peak  power  constraints  on 
individual  Trf\k)  can  also  be  introduced  as  in  [7]. 


III.  The  PIWF  Algorithm 


We  now  follow  the  same  principle  in  [7]  to  formulate  a 
cooperative  iterative  water  filling  scheduling  algorithm  for 
computing  'Kf\k). 

First  we  define  a  cooperative  utility  (or  objective)  function 
for  link  i  as  follows: 

K  K  Nt 

Ui  =  — 

k=l  k=l  1=1 

where  oji  can  be  any  useful  weight  for  link  i  (which  can  be 
a  product  of  the  distance  of  link  i,  a  penalty  factor  between 
zero  and  one  for  practical  coding  and  modulation,  and  other 
parameters),  and  Sf‘\k)  is  a  (non-negative)  price  factor  for 
channel  fc,  link  i  and  sub-stream  1.  The  first  term  in  is  the 
capacity  of  link  i,  and  the  second  term  in  Ui  is  a  penalty  term 
for  link  i.  It  is  the  second  term  that  provides  a  cooperation 
between  link  i  and  other  links.  Clearly,  the  rate  of  link  i  is 
directly  affected  by  the  power  distribution  ■nf\k)  for  all  I 
and  k.  At  the  same  time,  TTf\k)  for  all  I  and  k  also  cause 
interferences  to  other  links.  The  second  term  reflects  a  price 
that  inhibits  link  i  from  increasing  its  own  rate  selfishly. 
Note  that  each  penalty  component  to  the  utility  function  Ui 
is  purposely  chosen  to  be  simple,  which  is  the  product  of  a 
price  factor  and  the  corresponding  power  component.  If  Sl''\k) 
for  all  {k,l)  are  given,  and  link  i  wishes  to  find  7rP(fc)  to 
maximize  iti,  then  the  optimal  Trl''\k)  for  all  {k,l)  can  be 
shown  (see  the  appendix)  to  be 


where  [a;]+  =  max(a;,  0),  and  Pi  >  0  is  such  that 

^  (i)  _ 

^  ^  tt)  (k)  <  Kpi.  Furthermore,  if  there  is  no  positive 

k=l  1  =  1 

K  Nt  .  . 

Pi  such  that  (^)  =  Kpi,  then  choose  Pi  =  0.  The 

k=l  1=1 

above  is  what  is  called  priced-based  water  filling  [7].  If  there 
is  a  peak  power  constraint  on  7r|^^(fc),  the  optimal  7r|*^(fc)  is 
also  given  by  (7)  except  for  a  clipping  at  the  peak  value  [7]. 

In  order  to  determine  a  right  set  of  the  price  factors  6f\k) 
for  all  channels,  links  and  sub-streams,  we  set 


5f\k)  =  - 


Namely,  each  price  factor  for  link  i  is  set  to  be  equal  to 
the  negative  of  the  derivative  of  the  sum  of  the  weighted 
rates  for  all  other  links  with  respect  to  the  corresponding 
power  component.  This  rule  for  selecting  the  price  factors  is 
heuristically  sound.  It  was  also  obtained  in  [7]  by  comparing 
the  KKT  conditions  of  Mj  and  the  KKT 

conditions  of  max^O)  Ui  for  each  i. 

It  should  be  noted  that  the  domain  for  the  summation  in 


(8)  can  be  reduced  to  a  set  only  comprising  the  “neighboring” 
links  of  link  i.  A  link  j  is  considered  to  be  a  neighbor  of  link 
i  if  -^^-^jjj^u)jrj{k)  is  not  negligible,  or  equivalently  if  the 
channel  gain  between  the  transmitter  of  link  i  and  the  receiver 
of  link  j  is  not  negligible. 

We  next  provide  more  explicit  expressions  for  computing 
the  price  factors.  It  follows  from  (8)  that 


Sf\k) 

=  --^,T.^Aog,\i+p,{k)Gm 

dK  (k) 

=  {I  +  Pj{k)Gj{k)y^  Pj{k) 


=  -Hfyk)M-\k)  (  )  M-\k)H,yk) 

\dTi\pk)  J 

=  -HZ(k)M-\k)HtAk)(^^J^^PPk)^ 

■H^yk)M-\k)H,yk) 

=  a),  diag[nl^\k)  ■  •  • 

xoTil  [k)  / 

Nt 

=  -Tl^yk)diag  0  •  •  •  0  T  0  •  •  •  0  (10) 
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where  Tij{k)  =  Qf  {k)Hl^j{k)M~^{k)Hjj{k).  Since 
Tij{k)  is  proportional  to  Hij,  the  summation  in  (9)  needs 
to  be  over  the  set  of  all  j  where  j  ^  i  and  Hi  j  ^  0. 

The  computation  of  the  price  factors  6f  \k)  and  the  power 
components  T:f\k)  can  be  implemented  in  a  sequential  or 
concurrent  fashion  [7].  In  our  simulation,  we  will  consider 
the  sequential  implementation  as  follows. 

Step  1:  The  network  conducts  a  full  range  of  channel 
estimation.  By  the  end  of  this  process,  link  i  knows  Hj  j, 
Hji  and  Hi  j  for  all  j.  Also,  Pi{k)  =  0,  M~^{k)  = 
and  Gi{k)  =  are  set  for  all  {i,k),  and  a  counter 

is  initialized  at  f  =  0. 

Step  2:  For  each  new  value  of  t,  set  z  =  0. 

Step  3:  For  each  new  value  of  i,  link  i  makes  sure  that  it 
has  collected  the  most  recent  Pj{k),  M~^{k)  and  Gj{k)  for 
all  j  i  and  all  k  from  other  links.  Then,  link  i  computes  (or 
updates)  M~^{k)  and  Gi{k)  according  to  (4)  and  (3).  From 
Gi{k),  link  i  obtains  its  eigenvalues  xf\k)  for  all  {l,k)  and 
its  eigenvector  matrix  Qi{k)  for  all  k.  Then,  link  i  computes 
the  price  factors  6f\k)  for  all  {l,k)  according  to  (9)-(10). 
Finally,  link  i  updates  its  power  components  'Kf\k)  for  all 
((,  k)  according  to  (7)  and  consequently  updates  Pi{k)  for  all 
k. 

Step  4:  Set  i  =  i  +  1  and  go  to  Step  3  until  i  =  N. 

Step  5:  Set  t  =  t  +  I  and  go  to  Step  2  until  convergence. 

It  should  be  noted  that  in  practice  there  are  various  factors 
that  can  reduce  the  actual  network  throughput,  which  include 
non-ideal  coding  and  modulation,  channel  estimation  errors, 
quantization  errors,  etc.  However,  these  errors  can  be  con¬ 
trolled  by  the  network  designers  for  certain  applications.  The 
channel  estimation  issue  studied  in  [13]  is  among  many  useful 
considerations  in  practice.  There  are  also  many  different  ways 
for  channel  estimation.  The  PIWF  algorithm  clearly  consumes 
an  additional  amount  of  the  network  resources.  This  overhead 
should  be  relatively  small  compared  to  the  resources  used  for 
data  transmissions  (apart  from  the  information  exchanges  for 
in-network  computations  governed  by  the  PIWF  algorithm), 
or  otherwise  the  PIWF  algorithm  should  not  be  used.  How 
to  minimize  the  overhead  of  the  PIWF  algorithm  without 
reducing  its  performance  is  a  useful  topic  for  future  research. 
How  to  quantify  the  impact  of  channel  estimation  errors  and 
other  sources  of  errors  on  the  performance  of  the  PIWF 
algorithm  is  also  an  important  future  topic. 

IV.  Performance  Evaluation 
A.  Comparison  with  the  OPS  and  DI  Algorithms 

We  first  demonstrate  a  comparison  of  the  PIWF  algorithm 
with  the  OPS  algorithm  [4]  and  the  DI  algorithm  [1].  The  DI 
algorithm  is  simply  the  PIWF  algorithm  with  all  price  factors 
set  to  zero.  The  OPS  algorithm  is  a  centralized  algorithm 
which  maximizes  directly  the  network  throughput  using  a 
gradient-based  search.  We  will  choose  the  same  parameters 
as  used  in  [4],  i.e.,  the  signal  model  for  link  z  (z  =  1,  2, ...,  N) 
and  channel  k  {k  =  1, 2, ...,  K)  is 

N 

yi{k)  =  HijXi{k)+^  ^  HjjXj{k)  +  Wi{k)  (11) 


where  the  channel  matrices  Hj  i  G  (jN^xNt  ^j.g  assumed 
to  be  invariant  to  the  channel  index  k  and  consist  of  i.i.d. 
complex  Gaussian  random  entries  with  zero  mean  and  unit 
variance,  and  Wi{k)  is  complex  Gaussian  white  noise  vector 
with  zero  mean  and  the  covariance  matrix  NqI,  and  7  is  a 
parameter  used  to  model  the  total  interference  for  each  link. 
We  further  assume  that  N  =  6,  K  =  6,  Nt  =  Nr  =  2,  and 
lOlogiol^  =  20  where  p,  =  j^J2k=i  =  1- 

A  nominal  interference  to  noise  ratio  (INR)  in  dB  is  defined 
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as  lOlogj^g  Fig.  1  shows  the  sum  capacity  of  all  links 


Fig.  1.  Sum  capacity  in  bits/s/Hz  versus  nominal  INR  based  on  the  same 
data  as  in  [4]. 

averaged  over  R  =  100  channel  realizations  versus  the  INR 
in  dB  for  each  of  the  three  algorithms:  OPS,  PIWF  and  DI, 
all  with  the  unit  weight  oji  =  1.  The  averaged  sum  capacity  is 
computed  by  Y.r=i  Ek=i  Eti  where  rY\k)  is 

the  value  of  ri{k)  after  the  rth  run.  We  see  from  this  figure 
that  when  the  INR  is  small,  all  three  algorithms  yield  the 
same  capacity,  which  is  expected.  But  when  the  INR  is  large, 
the  capacity  differences  of  the  three  algorithms  become  large. 
The  OPS  has  the  largest  capacity,  the  DI  has  the  smallest 
capacity,  and  the  PIWF  is  somewhere  in  between.  This  figure 
also  shows  that  the  capacity  gap  between  the  OPS  and  the 
PIWF  is  still  quite  large  when  the  INR  is  large.  However,  the 
advantage  of  the  PIWF  over  the  OPS  is  that  the  computational 
cost  of  the  PIWF  is  much  smaller.  Indeed,  for  a  medium-sized 
network  as  shown  next,  the  complexity  of  the  OPS  algorithm 
became  so  overwhelming  that  its  simulation  was  not  feasible. 
Compared  to  the  DI,  however,  the  PIWF  is  more  costly  due 
to  the  computation  of  the  price  factors. 

B.  Throughput  of  the  PIWF  algorithm  for  a  larger  network 

We  now  consider  a  larger  network  of  MIMO  links  as 
illustrated  in  Fig.  2,  where  there  are  64  nodes  distributed 
randomly  (with  a  uniform  statistical  distribution)  in  a  square 
area,  and  32  links  {N  =  32)  are  formed  randomly  among 
these  nodes.  The  node  density  (i.e.,  the  number  of  nodes  in 
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Fig.  2.  A  network  of  64  randomly  distributed  nodes  and  32  links  for  random 
source-destination  pairs.  Node  density  is  one. 


Fig.  4.  A  network  of  64  randomly  distributed  nodes  and  32  links  between 
(approximately)  nearby  nodes. 


Fig.  3.  The  network  throughput  in  bits-meter/s/Hz/node  versus  the  traffic 
load  probability  where  K  =  1.  The  network  topology  is  shown  in  Fig.  2. 


a  unit  area)  is  one.  In  our  simulation,  we  had  100  network- 
wise  channel  realizations,  and  for  each  channel  realization  we 
formed  a  new  random  set  of  64  nodes  and  a  new  random 
set  of  32  links.  The  signal  model  for  this  network  is  given 
by  (1),  where  we  assume  that  Hij{k)  =  dij  is 

the  distance  between  the  transmitter  of  link  i  and  the  receiver 
of  link  j,  a  =  3  is  the  path  loss  factor,  and  Hij{k)  for  all 
{i,j,k)  consist  of  i.i.d.  complex  Gaussian  random  entries  of 
zero  mean  and  unit  variance.  The  noise  variance  Nq  is  such 
that  lOlogiQ  =  40  where  pi  =  1. 

Fig.  3  shows  the  network  throughput  in  bits- 
meter/s/Hz/node  versus  the  traffic  load  probability 
Here,  we  define  (  as  the  probability  that  each  link  (or 
link  i)  has  a  packet  to  transmit.  For  a  realization,  if  a 
link  has  no  packet  to  transmit,  it  will  not  take  part  in  the 
power  scheduling.  The  throughput  in  bits-meter/s/Hz/node  is 


computed  by 
C„ 


1 


R  N  K 


2NKR 


(12) 


r—1  i—1  k—1 


where  tOi  =  di^i.  For  this  figure,  only  a  single  channel  is  used, 

1. e.,  K  =  1.  We  see  that  the  capacity  increases  as  the  number 
of  antennas  increases,  which  is  expected.  Interestingly,  we  also 
see  that  for  I  <  Nt  =  Nr  <  4,  the  capacity  increases  almost 
linearly  with  the  number  of  antennas.  Another  useful  obser¬ 
vation  is  that  when  (  is  relatively  large,  the  capacity  of  the 
PIWF  for  such  a  network  topology  as  in  Fig.  2  decreases  as  ( 
increases.  This  observation  is  complementary  to  an  observation 
shown  in  [10]  that  the  shortest  distance  transmissions  are  the 
most  efficient  in  network  throughput  for  large  networks  unless 
the  traffic  load  is  low.  Indeed,  our  simulations  also  confirmed 
that  if  we  form  all  links  only  between  nearby  neighbors,  e.g., 
see  Fig.  4,  the  network  throughput  of  the  PIWF  algorithm 
becomes  an  increasing  function  of  (  for  the  entire  range 
0  <  c  <  1.  However,  a  disadvantage  of  the  nearby  neighbor 
transmissions  is  the  requirement  of  routing  and  flow  control 
protocols  between  sources  and  destinations,  which  adds  more 
complexity. 

Alternatively,  if  we  increase  the  number  of  bands,  as  shown 
in  Fig.  5  where  Nt  =  Nr  =  1,  the  capacity  at  higher  traffic 
load  can  be  improved. 

In  Table  I,  we  show  a  range  of  values  of  the  network 
throughput  under  the  PIWF  algorithm  for  1  <  Nt  =  Nr  <  6 
and  K  =  1,  2, 4, 8, 16, 32.  We  see  that  the  throughput  in  bits- 
meter/s/Hz/node  is  an  increasing  function  of  Nt  =  Nr  and 
K. 

It  is  useful  to  note  that  the  values  of  the  network  throughput 
in  bits-meter/s/Hz/node  (under  unit  node  density)  shown  in  Fig 
3,  Fig  5  and  Table  I  (especially  for  large  Nt  =  Nr  and  large 
K)  exceed  those  shown  in  [9]  and  [10].  However,  the  network 
considered  in  [9]  and  [10]  is  much  larger  than  the  one  in  Fig. 

2,  and  the  PIWF  algorithm  is  much  more  complex  than  those 


TABLE  I 

The  network  throughput  in  bits-meter/s/Hz/node  of  the  PIWF  algorithm  versus  the  number  oe  antennas  and  the  number  of 
CHANNELS.  N  =  Nt  =  Nr.  C  =  0.5.  THE  NETWORK  TOPOLOGY  IS  SHOWN  IN  FlO.  2. 
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Cmean 

K  =  1 

K  =  2 

K  =  A 

K  =  S 

K  =  16 

K  =  32 

N  =  1 

0.0513 

0.2096 

0.4005 

0.5553 

0.5897 

0.6114 

N  =  2 

0.2641 

0.5544 

0.9545 

1.1616 

1.1713 

1.1793 

N  =  3 

0.5149 

1.0609 

1.5562 

1.7353 

1.7607 

1.7682 

N  =  4 

0.7988 

1.5954 

2.0954 

2.2455 

2.2996 

2.3509 

N  =  5 

1.0952 

2.1062 

2.6106 

2.7446 

2.8000 

2.8100 

N  =  6 

I.4I40 

2.4604 

3.1140 

3.2915 

3.4740 

3.6033 

Fig.  5.  The  network  throughput  in  bits-meter/s/Hz/node  versus  the  traffic 
load  probability  where  Nt  =  Nr  =  1.  The  network  topology  is  shown  in 
Fig.  2. 


in  [9]  and  [10].  The  PIWF  algorithm  requires  and  exploits  the 
knowledge  of  all  interfering  sources  for  each  link.  For  the  case 
of  32  links,  4  bands  and  4  antennas  on  each  node,  for  example, 
it  took  about  5  minutes  (on  a  desktop)  for  the  PIWF  algorithm 
to  complete  the  computation  for  a  single  channel  realization 
under  full  load  condition.  For  the  network  size  (over  200 
nodes)  considered  in  [9]  and  [10],  the  simulation  of  the  PIWF 
algorithm  would  be  infeasible.  Although  the  neighborhood  of 
each  link  as  used  in  (8)  may  be  much  smaller  than  the  network 
size,  the  coupling  between  many  neighborhoods  makes  it 
difficult  for  the  PIWF  algorithm  to  converge.  In  fact,  we  also 
observed  that  if  the  radius  of  support  (with  respect  to  link 
i)  for  the  summation  in  (9)  reduces,  the  network  throughput 
produced  by  the  PIWF  algorithm  also  reduces. 

V.  Conclusion 

In  this  letter,  we  have  presented  an  extension  and  evaluation 
of  the  PIWF  algorithm  by  Wang,  Krunz,  and  Cui  [7]  for  an 
ad  hoc  wireless  network  of  MIMO  links.  The  extension  of  the 
PIWF  algorithm  to  the  MIMO  links  requires  each  link  to  have 
a  full  knowledge  of  the  channel  matrices  with  respect  to  all 
other  links,  and  the  computations  need  to  be  further  increased. 
The  study  shown  in  this  letter  provides  a  useful  perspective 
of  the  PIWF  algorithm  along  with  other  related  algorithms, 
which  is  not  available  in  [7].  This  work  also  suggests  that  there 
is  a  wide  range  of  complexity-performance  tradeoffs  in  the 


power  scheduling  algorithms.  Different  distributed  algorithms 
may  have  different  complexities  and  require  different  levels 
of  cooperations  between  links.  The  PIWF  algorithm  is  a 
distributed  algorithm  which  needs  a  high  level  cooperation 
between  links,  which  may  or  may  not  be  feasible  depending 
on  specific  applications  or  network  environment.  The  math¬ 
ematical  analysis  of  the  convergence  behaviors  of  the  PIWF 
algorithm  remains  an  open  problem. 


Appendix  A 
Proof  of  (7) 

The  optimization  problem  here  is  that  for  each 
i,  min^^U)(j,)j  — subject  to  —TTf\k)  <  0  and 

—  ^Pi-  'The  corresponding  Lagrangian 

function  is 


K  Nt 


fc=i 1=1 
K  Nt 


(13) 


vfc=l 1=1 


where  Ui  is  given  in  (6),  Q;f^(fc)  >  0  and  >  0  are  Lagrange 
multipliers.  The  K.K.T.  conditions  [11]  for  the  optimal  Tif\k) 
are  given  as  follows: 
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dnf\k) 
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dnf\k) 


K  Nt 

EE 

k^l 1^1 
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(15) 
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(16) 
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We  also  know 


— II) — 

d4\k) 

-4— log2|/  +  ^*(fc)G,(fc)| 
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=  (logs  e)  - 


Then,  using  (20)  in  (19)  yields 


-  Ti)-  (21) 
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where  [x]’*'  =  max(a;,  0)  is  chosen  to  accommodate  those 
{l,k)  where  nfl^k)  =  0.  Namely,  for  those  {l,k)  where 
Trf{k)  =  0,  there  must  be  df\k)  >  0  such  that  (14)  holds, 
which  is  equivalent  to  the  operation  [x\^  used  in  (21).  From 
(17)  and  (21),  the  parameter  /3*  >  0  should  be  such  that 


10g2  e  ^  1 

hh[pi  +  4‘\k)  xf{k) 


<  Kp^ 


Note  that  because  of  the  price  factors  (5)  (fc),  it  is  possible 
that  the  equality  in  (22)  may  not  need  to  be  achieved  for  some 
cases.  Indeed,  since  Trf{k)  as  in  (21)  is  a  decreasing  function 


K  Nt 

of  (3i,  if  E  E 


uji  log2  e 


Pi  is  zero  according  to  (18). 


<  Kpi,  then  the  optimal 


